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We present results for the ordered sequence of first passage times of arrival of TV random walkers 
at a boundary in Euclidean spaces of d dimensions. 



. 

It is customary to assume that the important behaviors of a group of independent events can all be characterized 
by studying the behavior of one such event. For instance, one focuses on the mean time that it takes a single random 
walker to arrive somewhere even when there are many such random walkers in the system provided the walkers are 
independent. However, it is clear that there are situations in which one might be interested, for example, in the mean 
time for the first of the walkers to arrive somewhere or, more generally, in the ordered sequence of a particular outcome. 
The statistics of ordered outcomes clearly depends on the number N of events even when these are independent. For 
example, the mean time for the first of /V walkers to arrive somewhere must clearly decrease with increasing N. The 
interest in so-called order statistics has grown with the development of experimental techniques that make it possible 
to follow single events on the microscopic scale. 

One way to characterize the behavior of N random walkers is to focus on the order statistics for arrival at a 
' O ■ boundary, and, in particular, the mean and variance for the first arrival of the first, second, third, etc. walker out of a 
set of N. Most of the early literature dealt with mean arrival times in one dimension More recent literature has 

extended these concepts to fractal lattices [pj-pt and has attempted to include higher-dimensional Euclidean lattices 
as well ||[7j] . Some of these recent efforts [Qhave relied on a parallel literature about another quantity of interest in 
these problems, namely, the distinct number of sites visited by N walkers [||. In this approach the first arrival time 
statistics are obtained through a conjectured relation between arrival times and distinct number of sites visited. Yer 
another approach |J relies on scaling arguments that are said to be independent of the underlying environment and 
can be applied in any dimension to both ordered and disordered structures. However, the results obtained in this work 
£f} \ do not agree with any of the other published results (nor with those obtained herein), even in the well-established 
i one-dimensional problem. This difference is apparent in the leading term of the result and also in the form of the 
series that follows the leading term. 

In addition to order statistics and distinct number of sites visited, the behavior of N random walkers has also been 
characterized in terms of the maximal excursion P,|lO|]. 

In order to complete the iV-independent-walker panorama, in this report we present results for the ordered sequence 

S of first passage times of arrival at a boundary in Euclidean spaces of d dimensions. This completion is made possible by 
a result reported in the maximal excursion literature [|| . The leading terms obtained here agree with those obtained 
■ earlier and therefore confirm the conjecture made in that work (at least to leading order). 

The calculation of ordered sequences of first passage time moments is rather elaborate but has been laid out in 
O ■ detail in a number of previous papers §-|j. These steps are easy to state but complex to carry out. The mth moment 
' of the jth passage time (that is, of the first arrival of the jth particle out of N) at a distance r from the origin is 



(*i»= / t m ^ N {r,t)dt. (1) 



Here ipj t N(r,t) is the probability density of the time it takes the jth particle to first reach the given distance r. 
This quantity can in turn be expressed in terms of the first-passage time density to this distance of a single particle, 
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^(r,*) = ^i,i(r,f) fcUl 



V^(r,t) = N\/[(N-j)l(j - l^M^-V^Ml - h(r,t)]^- J ' 



(2) 



Here h(r,t) = ^j(r,T)dr is the mortality function, i. e., the probability that a single diffusing particle has reached 
the distance r during the time interval (0,i). Thus, knowledge of the function ^){r,t) ( or h(r,t)) in principle allows 
the full evaluation of first passage time moments - provided one can carry out the integral in Eq. (|]J). Alternatively, 
one can calculate the generating function of the moments, 
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(3) 



and from this obtain the moments via a Taylor series expansion in z. 

Herein, of course, lies the difficulty of the problem: neither can the integral in Eq. (Q) be done easily, nor can 
the generating function be calculated easily. Indeed, the literature is based on rather elaborate expansions and/or 
renormalization group procedures to obtain the generating function, methods that rely on there being a very large 
number of walkers (N large) and on the walker of interest being one of the first few (j <C TV) so that short times 
dominate the moments. The formal results are similar from one system to another (Euclidean, fractal), the differences 
arising because of the differences in the mortality function and, specifically, in the short time behavior of this function. 

The mortality function (and particularly its short time behavior) for one-dimensional systems was calculated in the 
early work on order statistics [^|J|] . For certain deterministic fractal geometries this function was considered in more 
recent work [0,0. In all of these, the short-time behavior of h(r, t) has a "universal" form. It depends only on the 
combination t/r 2 (which we denote simply as t, understanding that it is the distance-scaled time) and is given bv 



(4) 



where the tilde just stresses the fact that this is an asymptotic result for t — ► 0. The constants A, fi, /3, 7, and S vary 
from one system to another. 

Although mortality functions for Euclidean systems have long been known and used in the random walk literature 
for a variety of problems [p| JT^Jl^| , its short-time behavior for systems of dimension d > 1 has not been incorporated 
into the order statistics context. Perhaps not surprisingly, it turns out to fit the pattern Eq. (Q) and therefore the 
existing theories can directly be applied these systems. In particular, the time Laplace transform of the mortality 
function h(r, s) — Ch(r,t) (indicated by the same symbol as the function but with the argument t replaced with the 
Laplace variable s) is 



sh(r, s) 



2 l-d/2 (2dr 2 s) 



where I v is a modified Bessel function of order u, and 



!d/2-i{z) 



(d- l)(d-3) 
8~z 



+ 0(z- 2 ) 



(5) 



(6) 



for z ^> 1. In writing (0) an implicit choice of the diffusion coefficient D in the standard mean squared displacement 
relation (r 2 ) = 2dDt has been made, namely, D = l/2d, so that (r 2 ) = t in each Euclidean dimension. The inverse 
Laplace transform of s v exp(— as 1 / 2 ) is 



-1 / v ~as 1 ' 2 \ 1 -a 2 /(8«) n f G \ 



(7) 



where D n {x) is a parabolic cylinder function (or Whitaker's function). Asymptotic expansion of this function for 
x \n\ |L5| leads for r 2 ft 3> 1 to 



h(r, t) 



dr 2 



r(d/2) V 2t 
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(8) 



or, in terms of the scaled time t/r 2 — > t, to the small t result 



h{t) = At»e- d /( 2t \l + hit) 



(9) 
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with 
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This is precisely of the form (|4|). 

Since the steps leading to the moment expansions are well documented in the literature once the form (Q) has been 
established [|f-^|, we only present the results. For the m-th moment of the first passage time of the first of N ^> 1 
particles we find (the distance r does not appear explicitly because all times are scaled) 
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where 



(13) 



An = 



r(d/2) 

The m-th moment of the first passage time of the j-ih particle with j N is 

jm J — 1 



(Cat) » <*5» + 
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where 



A„(m) = 1 + 
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(14) 
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and Si(n) is a Stirling number of the second kind |l5|. In particular, the variance a 2 N = (t 2 N ) — (tj.N) 2 can be 
obtained from Eq. (pi) with (pi): 
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To check on the range of validity of these results it would be desirable to carry out direct computer simulations 
involving a very large number of walkers, an exercise that is costly. Another approach is to integrate numerically 
the exact Eq. ([!]) with Eq. (J2j> and the exact mortality function in the integrand for values of t beyond the range of 
validity of the short-time expansions. The exact survival probability in one dimension is fl2]| 
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(18) 



In two dimensions fl2|| 



n~^0 X 0nJl( x 0n) 



exp 



4r 2 



(19) 
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where J\ is the first Bessel function and xo n is the n-th root of 1 
Hollingsworth distribution Jl3| 

S(r,t) = l- — 2 ex P 

^ ^ n— — 00 

To carry out the integration it is convenient to integrate by parts and so that 

/VI f°° 

(t™ N ) = (tftr) + »' ( V j!!y! / * m-1 [l - S(r, t)]'S(r, t) N ^dt (21) 

where 

/•OO 

(t? N )=m t m - l S{r,t) N dt. (22) 
./o 

Note that these integrations become rather awkward for large values of N. 

The numerical results and comparisons with asymptotic expressions are presented graphically in a set of figures. In 

Fig. [| we present scaled results for the N dependence of (^iV) 1 ^™, that is, the 1/mth power of the mth moment of 
the first passage time of the first of N particles. The exact results obtained from the integration ( p2|) are indicated 
by symbols and the theoretical results of various orders by lines. Circles denote exact results for d — 1, triangles for 
d = 2, and squares for d = 3. The first panel shows results for m = 1, that is, for the mean first passage time of the 
first particle to the desired boundary; the second panel presents the second moment (m = 2), and the third panel the 
third moment. The moments on the ordinates are scaled so that the scaled moment approaches 1 as 1/ In N — » 0. The 
integrations were performed for N = 2 3 , 2 4 , • • • , 2 30 . The dotted curves correspond to the asymptotic results Eq. ( p"3| ) 
to zeroth order, that in, only the first term in the series. The dashed lines include two terms, the solid lines three. 
Clearly the convergence to the exact results improves with order retained, but more slowly with increasing dimension. 
In any case, the deviations even at order 2 are clear on the scale of the figure at around N = 100. 




1/lnN l/\nN MnN 

FIG. 1. Scaled moments for first arrival of the first of N particles at a prescribed boundary as a function of 1/lnN. First 
panel: first moment (i. e., mean first passage time); second panel: second moment; third panel: third moment. Numerical 
results are indicated by circles (d = 1), triangles (d = 2), and squares (d = 3). Asymptotic results (cf. Eq. (|l3|)) to zeroth 
order, first order, and second order are shown by dotted curves, dashed curves, and solid curves respectively. 

The two panels shown in Fig. [2] deal with the variance o\ N = {t\ N ) — (ti,jv) 2 - The leading term in the large-TV 
analysis of the variance is given in Eq. (^) ; for the variance we have only calculated this leading (zeroth order) term. 
In the first panel we present the 1/ In A^-dependent behavior of the variance scaled in such a way that the leading 
term in the theoretical expression gives unity for all the dimensions d = 1,2,3 (indicated by the dotted line in the 
figure). The circles (d — 1), triangles [d = 2), and squares (d = 3) are the numerical results obtained from explicit 
integration. A number of interesting observations are apparent from this panel. First, we see that the leading term 
in the expansion of the variance leads to adequate results for some range of N only in one dimension. In two and 
three dimensions even for extremely large it is necessary to go beyond zeroth order to obtain adequate analytic 
results. It is thus clear that the correction terms for the variance are much more important than for the first passage 
time moments. This can be seen from our asymptotic expansions since the ratio of first to zeroth order terms goes 
as ln(lniV)/liiiV for the moments but as [In (In N)} 3 / In N for the variance. In fact, the second actually grows with 



tie Bessel function Jq. For d — 3 one obtains the 
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increasing N at first, becoming larger than unity, and only begins to decrease for extremely large values of N (of order 
2 30 ~ 10 9 ). This in turn leads to the "anomaly" observed in the first panel and enlarged in the inset, namely, that 
the exact variance actually crosses and becomes larger than the zeroth order theoretical one before settling down to 
the asymptotic value at extremely large N. (We only carried the numerical calculations to these very large N values 
for one-dimensional systems and in fact reach the limit of numerical reliability in that region.) 

The second panel in Fig. \3 shows the same information as the first but plotted in a different way to stress other 

— 1/2 

features. It is simply a plot of the a l ^ vs In Y comparing numerical (symbols) and zeroth order asymptotic (lines) 
results in one, two, and three dimensions. Again, it is clear that for d = 1 the zeroth order asymptotic result is quite 
good but for higher dimensions it is not adequate. 
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FIG. 2. N dependence of the variance for first arrival of first of N walkers at a prescribed boundary. Symbols and lines are 
as in Fig. |l[ First panel: variance scaled so that the zeroth order asymptotic result is unity in all dimensions (indicated by the 
dotted line). Inset: detail of the d — 1 numerical results. Second panel: same results, unsealed. 



In conclusion, we have filled some missing pieces in the mosaic of results for the order statistics of N independent 
random walkers in d-dimensional lattices. These results complement and confirm previous conjectures, indicate that 
large- TV asymptotic expansions converge more rapidly in lower dimensions, and that convergence of higher cumulants 
such as variances is more problematic than that of first passage time moments. 
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